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^ : 1 Introduction 

o. 

Hyperbolic conservation laws posed on manifolds arise in many applica- 
tions to geophysical flows and general relativity Recent work by the author 
1^ I and his collaborators attempts to set the foundations for a study of weak 

solutions defined on Riemannian or Lorentzian manifolds and includes an 



investigation of the existence and qualitative behavior of solutions. The 
^ . metric on the manifold may either be fixed (shallow water equations on 

F^ ', the sphere, for instance) or be one of the unknowns of the theory (Einstein- 

Euler equations of general relativity). This work is especially concerned 
with solutions and manifolds with limited regularity. We review here re- 
sults on three themes: (1) Shock wave theory for hyperbolic conservation 
laws on manifolds, developed jointly with M. Ben-Artzi (Jerusalem); (2) 
Existence of matter Gowdy-type spacetimes with bounded variation, de- 
veloped jointly with J. Stewart (Cambridge). (3) Injectivity radius estimates 
for Lorentzian manifolds under curvature bounds, developed jointly with 
B.-L. Chen (Guang-Zhou). 



2 Conservation laws on a Riemannian manifold 

In the present section, (M", g) is a compact, oriented, n-dimensional Rie- 
mannian manifold. As usual, the tangent space at a point x G M" is 
denoted by TM" and the tangent bundle by TM" := [JxeU" TM", while 
the cotangent bundle is denoted by T*M" = r*M". The metric structure is 
determined by a positive-definite, 2-covariant tensor field g. 

A flux on the manifold M" is a vector field / = fx{u) depending smoothly 
upon the parameter u. The conservation law associated with / reads 

dtu + div^ {f{u)) = 0, (2.1) 

where the unknown u = u{t,x) is defined for t > and x G M" and the 
divergence operator is applied to the vector field x i-> fx{u{t,x)) G T^M". 
We say that the flux is geometry-compatible if 

div gf^iu) = 0, w G R, X G M". (2.2) 

We propose to single out this class of conservation laws as an important 
case of interest, which leads to robust U estimates that do not depend on 
the geometry of the manifold. The equations arising in continuum physics 
do satisfy this condition. 

Equation (|2.1|) is a geometric partial differential equation which de- 
pends on the geometry of the manifold only. All estimates must take a 
coordinate-independent form; in the proofs however, it is often convenient 
to introduce particular coordinate charts. We are interested in solutions 
u G L"'(R+ X M") assuming a prescribed initial condition Uq g L°°{M."): 

u{0,x) = Uo{x), xeM". (2.3) 

We extend Kruzkov theory of the Euclidian space IR" (see [8J) to the case 
of a Riemannian manifold, as follows. 

Let / = fx{u) be a geometry-compatible flux on the Riemannian man- 
ifold (M,g). A convex entropy/entropy-flux pair is a pair {U,F) where 
LT : IR — > R is convex and F = Fx(Q) is the vector field defined by 



(w) := I du' 



F:,{u) := du'U{u')du'fx{u')du' , w G R, x G M" 



Given Uq e L°°(M") a function u e L°°^+,L°°(M")) is called an entropy 
solution to the initial value problem (|2.1[) , (|2.3|) if the following entropy 
inequalities hold 



J Jr 



(Lr(u(f, x)) <9f 0(f, x) + g.(f x(w(f, x)), grad^ 0(f, x))) rfy^(x)df 



R+xM" 

+ I L[(Mo(x))0(O,x)dy^(x)>O, 
Jm« 

for every convex entropy/entropy flux pair {U, F) and all smooth function 
6 = d{t,x) > compactly supported in [0, oo) x M". 

Theorem 2.1 (Well-posedness theory on a Riemannian manifold. I). Let 

f = fx{u) be a geometry-compatible flux on a Riemannian manifold (M",g). 
Given Uq e L°°(M") f^ere exists a unique entropy solution u G L'^(R+ x M") to 
the problem (I2.1P - (|2.3P . Moreover, for each 1 <p <oo, 

\\^{i)\\LP(M";dVg) ^ \\Uo\\LP{M";dVg), t ^ 1R+/ 

and, given two entropy solutions u, v associated with initial data Uq, Vq, 

\\v{t) - M(OllLi(M";rfy^) ^ Ibo - WollLi(M«;rfyj,)' ^ ^ ^+- 

The framework proposed here allows us to construct entropy solutions 
on a Riemannian manifold via the vanishing diffusion method or the finite 
volume method IHIll. Following DiPema [6J we can introduce the (larger) 
class of entropy measure-valued solutions {t, x) G R+ x M" i— > Vt x- 

Theorem 2.2 (Well-posedness theory on a Riemannian manifold. II). Let 

f = fx{u) be a geometry-compatible flux on a Riemannian manifold {M.",g). Let 
V be an entropy measure-valued solution to (l2.1P - (l2.3P /or some uq g L'"(M"). 
Then, for almost every {t,x), Vt,x = 5u(t,x)r where u G L"'(R+ x M") is the unique 
entropy solution to the problem. 

Finally we can relax the geometry compatibility condition and consider 
a general conservation law associated with an arbitrary flux /. More 
general conservation laws solely enjoy the L^ contraction property and 
leads to a unique contractive semi-group of entropy solutions. 



Theorem 2.3 (Well-posedness theory on a Riemannian manifold. III). Let 

f = fx{U) be an arbitrary (not necessarily divergence-free) flux on {M.",g), satis- 
fying the linear growth condition 

\fx{u)\g < 1 + \u\, i/ G R, X G M". 

Then there exists a unique contractive, semi-group of entropy solutions Uq g 
L^{M") H^ u{t) := StUo G L^(M") to the initial value problem (IJTID , (|2J] ). 

For the proofs we refer to fT]. See fT] for the convergence of the finite 
volume schemes on a manifold. Earlier material can be found in Panov [12l 
(n-dimensional manifold) and in LeFloch and Nedelec [lOJ (Lax formula 
for general metrics including the case of spherical symmetry). 



3 Conservation laws on a Lorentzian manifold 

Motivated by the application to general relativity, we can extend the the- 
ory to a Lorentzian manifold. Let (M"'^^,g) be a time-oriented, (n -I- 1)- 
dimensional Lorentzian manifold, g being a metric tensor with signature 
{—,+,..., +). Tangent vectors X can be separated into time-like vectors 
{g{X, X) < 0), null vectors (g(X, X) = 0), and space-like vectors {g{X, X) > 0). 
The null cone separates time-like vectors into future-oriented and past- 
oriented ones. Let V be the Levi-Cevita connection associated with the 
Lorentzian metric g. 

A flux on the manifold M""^^ is a vector field x h^ fx{u) G T^M"'^^, 
depending on a parameter u G R. The conservation law on (M""^^,^) 
associated with / is 

div^ (/(u)) =0, u: M"+^ -^ R. (3.1) 

It is said to be geometry compatible if 

div^/;,(fi) = 0, QgR, xgM"+\ (3.2) 

Furthermore, / is said to be a tim.e-like flux if gx[dufx{u),dufx{ii)) < 0, 
X G M"+\ Q G R. 

Note that our terminology here differs from the one in the Riemannian 
case, where the conservative variable was singled out. We are interested 



in the initial-value problem associated with (|3.1|) . We fix a space-like hy- 
persurface [Kq c M""^^ and a measurable and bounded function Uq defined 
on IKq. Then, we search for u = u{x) e L°°(M"'^^) satisfying (|3.1|) in the 
distributional sense and such that the (weak) trace of u on "Kq coincides 
with Uq: 

U\:ko = Uq. (3.3) 

It is natural to require that the vectors dufxiu) are time-like and future- 
oriented. 

We assume that the manifold M""^^ is globally hyperbolic, in the sense 
that there exists a foliation of M""^^ by space-like, compact, oriented hy- 
persurfaces !Kt {t G IR): M"+^ - UfeR ^f- ^^ly hypersurface 'Ktg is referred 
to as a Cauchy surface in M""^^, while the family "Kt {t G R) is called an 
admissible foliation associated with "Kt^. The future of the given hyper- 
surface will be denoted by M"+^ := [Jt>o "^t- Finally we denote by n^ the 
future-oriented, normal vector field to each "Kt, and by g'^ the induced met- 
ric. Finally, along "Kt, we denote by X^ the normal component of a vector 
field X, thus X*:=g(X,nO. 

A flux F = Fx{u) is called a convex entropy flux associated with the 
conservation law (|3.1|) if there exists a convex LT : R — » IR such that 

Fx{u)= I duU{u')dufx{u')du', xgM"+^mgR. 
Jo 

A measurable and bounded function u = u{x) is called an entropy solution 
of the geometry-compatible conservation law (|3.1 |) -(|3.2 |) if 



r g{F{u),grad gd)dVg+ f go{F{uo),no)d:^dVg^>0. 

for all convex entropy flux F = Fx{u) and all smooth > compactly 
supported in M""^^. 

Theorem 3.1 (Well-posedness theory on a Lorentzian manifold). Consider 



a geometry-compatible conservation law (I3.1P - (13.2P posed on a globally hyperbolic 
Lorentzian manifold M"'^^. Let ^o be a Cauchy surf ace in M."'^^, and Uq : CKq — > R 
be measurable and bounded. Then, the initial-value problem (I3.1P - (I3.3P admits a 
unique entropy solution u = u{x) G L°°{M."'^^). For every admissible foliation 'Kt, 
the trace Uo^^ exists and belong to L^CKt), and ||f'(w:K(llLi(jc,) is non-increasing in 



time, for any convex entropy flux F. Moreover, given any two entropy solutions 
u,v. 



||/(t/^,)-/(i;|:K,)llLM^,) 



is non-increasing in time. 



We emphasize that, in the Lorentzian case, no time-translation property 
is available in general, contrary to the Riemannian case. Hence, no time- 
regularity is implied by the L^ contraction property. 

4 Existence of matter Gowdy-type spacetimes 
with bounded variation 

Vacuum Gowdy spacetimes are inhomogeneous spacetimes admitting two 
commuting spatial Killing vector fields. The existence of vacuum space- 
times with Gowdy symmetry is well-known and the long-time asymptotics 
of solutions have been found to be particularly complex. In comparison, 
much less emphasis has been put on matter spacetimes. Recently, LeFloch, 
Stewart and collaborators HI HI]] initiated a rigorous mathematical treat- 
ment of the coupled Einstein-Euler system on Gowdy spacetimes. The 
unknowns of the theory are the density and velocity of the fluid together 
with the components of the metric tensor. The existence for the Cauchy 
problem in the class of solutions with (arbitrary large) bounded total varia- 
tion is proven by a generalization of the Glimm scheme. Our theory allows 
for the formation of shock waves in the fluid and singularities in the ge- 
ometry. The first results on shock waves and the Glimm scheme in special 
and general relativity are due to Smoller and Temple |]T4| (flat Minkowski 
spacetime) and Groah and Temple |]Zi (spherically symmetric spacetimes). 
The novelty in |]3] E]] is the generalization to a model allowing for both 
gravitational waves and shock waves. 

The metric is given in the polarized Gowdy symmetric form 

ds^ = ^ i-dt^ + dx^) + e^^ {e^' dy^ + e-^' dz^), (4.1) 

where the variables a,h,c depend on the time variable t and the space 
variable x, only. We consider Einstein field equations G^f^ = kT"^ for 
perfect fluids with energy density /j > and pressure p = [j.cl. Here, the 



sound speed Cs is a constant with < Cs < 1 and G"^ denotes the Einstein 
tensor, while k is a normalization constant. 

The 4-velocity vector u" of the fluid is time-like and is normalized to 
be of unit length and we define the scalar velocity v and relativistic factor 
^ = ^(z;)by(M") = e"" ^ (1, z;, 0, 0) and ^ = {l-v'^)~^^'^. The matter is described 
by the energy-momentum tensor T"f^ = {[i +p) u"u^ + p g"^, from which we 
extract the fields t, S and L: 

T°°= e-^{{^ + p) e - p) =:e-^T, 

After very tedious calculations we arrive at the constraint equations 

lutht + 2a^h^ + h] - 2h^^ - 3hl - c^^ - cl =k e^ %, 
-2 atbx-2a^bt + 2 btx + 2btbj + 2ctCx =k e^ S, 

and the evolution equations 

att -a^x =b]-bl-c] + cl + - e^" (-t + L - 2p), 

bn-bxx =-2b^ + 2bl + '^e^iT-Z), 
Ctt-Cxx = -2btCt + 2bxCx. 
The evolution equations for the fluid, V^T"^ = 0, are the Euler equations 

Tt + Sx = Ti, St + 'Lx = Ti, 

in which the source terms Ti, T2 are nonlinear in first-order derivatives of 
the metric and fluid variables. 

We propose to reformulate the Einstein-Euler equations in the form of 
a nonlinear hyperbolic system of balance laws with integral source-term, 

in the variables {ii,v) and w := (at,ax,(it,(ix,Ct,Cx), where (i = e^^. It is 
convenient to also set a = e^". The functions a, b (and a,(i) are determined 
by 

a{t, x) = e^^^'^\ a{t, x) = j_^ W2{t, y) dy, 
b{t, x) = ^ In j6(f , x), j6(f, x) = 1 -I- J_^ w^it, y) dy. 



Obviously, we are interested in solutions such that (i remains positive. 

The equations under consideration consist of three sets of two equa- 
tions associated with the propagation speeds ±1, the speed of light (after 
normalization). The principal part of the fluid equations are the standard 
relativistic fluid equations in a Minkowski background, with wave speeds 
A+ = {v± Cs)/(1 + V Cg). To formulate the initial-value problem it is natural to 
prescribe the values of [i, v, w on the initial hypersurf ace at f = 0, denoted 
by(fi°,z;°,z(;°). 

Our main result is: 

Theorem 4.1 (Existence of Gowdy spacetimes with compressible matter). 
Consider the {[j.,v,w)-formulation of the Einstein-Euler equations on a polarized 
Gowdy spacetime with plane-symmetry. Let the initial data {[i'^,v°,w'^) be of 
bounded total variation, TV{[i°,v°,w'^) < oo, satisfying the constraints, and 
suppose that the corresponding functions a^,b^ are measurable and bounded, 
sup |(tt°, fc°)| < oo. Then the Cauchy problem admits a weak solution [j., v, w such 
that for some increasing C{t) 

TV{^, V, w){t) + sup \{a, b){t, ■)\ < C{t), t > 0, 

and are defined up to a maximal time T^oo.IfT<oo then either the geometry 
variables a, b blow up: limf^^fsupj^ \a{t, ■)\ + \b{t, ■)\\ = cxj, or the energy density 
blows up : lirat^T supj^ \^{t,-)\ = oo. 

Hence, the solution exists until either a singularity occurs in the geom- 
etry (e.g. the area (i of the 2-dimensional space-like orbits of the symmetry 
group vanishes) or the matter collapses to a point. To our knowledge this 
is the first global existence result for the Euler-Einstein equations. 

If a shock wave forms in the fluid, then [j-,v will be discontinuous 
and, as a consequence, Wsx and Zf4x might also be discontinuous. In fact. 
Theorem 5 allows not only such discontinuities in second-order derivatives 
of the geometry components (i.e. at the level of the curvature), but also 
discontinuities in the first-order derivatives which propagate at the speed 
of light. The latter correspond to Dirac distributions in the curvature of 
the metric. 



5 Lower bounds on the inj ectivity radius of Loren- 
tzian manifolds 

Motivated by the application to spacetimes of general relativity and by 
earlier results established by Anderson [2J and Klainerman and Rodni- 
anski flSl, we investigate in [5] the geometry and regularity of (n + 1)- 
dimensional Lorentzian manifolds {M,g). Under curvature and volume 
bounds we establish new injectivity radius estimates which are valid ei- 
ther in arbitrary directions or in null cones. Our estimates are purely local 
and are formulated via the "reference" Riemannian metric ^j associated 
with an arbitrary future-oriented time-like vector field T. 

Our proofs are based on suitable generalizations of arguments from 
Riemannian geometry and rely on the observation that geodesies in the 
Euclidian and Minkowski spaces coincide, so that estimates for the ref- 
erence Riemannian metric can be carried over to the Lorentzian metric. 
Our estimates should be useful to investigate the qualitative behavior of 
spacetimes satisfying Einstein field equations. 

We state here one typical result from |5l encompassing a large class of 
Lorentzian manifolds. Fix a point p G M, and let us assume that a domain 
Q c M containing p is foliated by spacelike hypersurfaces Zf with normal 
T, say Q = IJfe[-i,i] ^f- Assume also that the geodesic ball B^gipA) c Zq 
is compactly contained in Lq. Consider the following assumptions where 
Kq, Ki, K2 and Vq are positive constants: 



Ko.-ll 



2 



8 



<l/Ko inQ, (5.1) 



\LTg\g, < Ki in n, (5.2) 

IRm^l^^ < ^2 in n, (5.3) 

volg{B^,{p,l))>Vo, (5.4) 



We prove in ||5l : 



Theorem 5.1 (Injectivity radius estimate for Lorentzian manifolds). Let 
(M, g) be a Lorentzian manifold satisfying (I5.1P - (I5.4P at a point p e M. Then, there 
exists a positive constant z'o depending on the foliation bounds Ko,Ki, curvature 
bound K2, volume bound Vq, and dimension n so that the injectivity radius at p is 
bounded below by z'o, that is Ini{M,g,p) > Zq. 
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